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Abstract. Recently, new data on antiprotons and positrons from PAMELA, e +e + spectra from 
ATIC, FERMI and HESS up to TeV energies all indicate deviations from expectations, which has 
caused an interesting mix of new explanations, ranging from background, standard astrophysical 
sources to signals from dark matter (DM) annihilation. Unfortunately, the excess in positrons is not 
matched with obvious excesses in antiprotons or gamma rays, so a new class of DM scenarios with 
"leptophilic" WIMP candidates have been invoked. On the other hand, the increase in the positron 

, , fraction, which could have had any spectral shape for new physics, matches well the shape expected 

from proton background. 
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P INTRODUCTION 

The quest for understanding the nature of the elusive dark matter has received a new 
boost by the rising positron fraction observed by the PAMELA space-borne experi- 

i—i ment QJ in combination with rather hard spectra of electrons and positrons above ex- 

pected background by the FERMI satellite O, the ATIC balloon experiment 01 and 
the H.E.S.S. earth-bound Cherenkov telescope [4J. These excesses were not accompa- 
nied by an obvious excess in diffuse gamma rays in the halo at mid-latitude Q nor in 
antiprotons [6], which has led to speculations about a new class of "leptophilic" dark 
matter candidates HIEL which fit the data BUl. On the other hand astrophysical explana- 
tions, like nearby pulsars IIISlIIIlIEllISQIIElIISIiaor nearby SNRs El II [13, 

O provide viable explanations. 

> In this short summary we first discuss the data and then the possible explanations in 

terms of astrophysical sources or dark matter annihilation. Given the explosion of papers 
it will not be possible to discuss all ideas. Recent reviews on dark matter can be found 
inRefs. lEHI and rf2TTl. 
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PAMELA DATA ON POSITRONS AND ANTIPROTONS 

PAMELA (Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics) 
is a satellite-borne experiment carried on board of a Russian earth observing satellite, 
which was launched in June, 2006. Its main components are i) a magnetic spectrometer 
consisting of a silicon tracker inside a 45 cm high permanent magnet with a square 
inner bore of 16x13 cm 2 and a field of 0.4T; ii) a 24x24 cm 2 Si-tungsten calorimeter 




with a thickness of 16.3 Xo; iii) a time-of-fiight counter iv) an anti-coincidence counter; 
v) a neutron counter. The main purpose of the experiment is the measurement of the 
antiproton and positron components of cosmic rays in a energy range of 50 MeV to a 
few 100 GeV with high statistics. 

The main surprise of the first two years of data taking was the observation of an 
increasing positron fraction, defined as the fraction e + /e + + e~, above 10 GeV, while 
standard propagation models, like GALPROP, predict a continuously decreasing fraction 
QT The data up to 100 GeV are shown in Fig.[T^ and agree with previous measurements 
up to 10 GeV by HEAT J221 and AMS-01 fl23]|7ln contrast, the antiprotons did not show 
any particular feature, as shown by the antiproton/proton ratio in Fig. O. 

ELECTRON PLUS POSITRONS FROM ATIC, FERMI AND HESS 

The Advanced Thin Ionization Calorimeter (ATIC) is a balloon-based experiment de- 
signed to study the spectra of cosmic ray protons, light nuclei, and electrons. It measures 
electromagnetic showers in a 16 Xo thick calorimeter. The spectrum of the summed 
fluxes of positrons and electrons start to become harder above 300 GeV with an abrupt 
decrease at around 600 GeV, as shown in Fig. [2^ [3J. The rather steep decline, observed 
in 3 independent flights, has led to speculations about new physics, especially the pos- 
sibility of the dark matter annihilation of Kaluza-Klein type WIMPs. These particles 
occur naturally in theories of universal extra dimensions as reviewed by Hooper and 
Profumo ll23Tl . They can decay into lepton pairs, which would lead to a peak in the spec- 
trum smeared with a low energy tail from final state radiation. Given that the excess stops 
around 800 GeV this would require WIMP masses in this range. However, the bump was 
not confirmed by data from the FERMI telescope 0, a gamma-ray telescope launched 
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FIGURE 2. Left: Spectrum of electrons+positrons from FERMI and ATIC. From Abdo et al. (2:| . Right: 
as left, but for different experiments and including new data from H.E.S.S. From Aharonian et al. |4). 



in June 2008 by the NASA. It consists of a silicon tracker including a converter to con- 
vert gamma rays into electron pairs, whose energy is measured in a CsI(TI) scintillating 
calorimeter with a thickness of 8.6 Xrj. The FERMI data showed a smooth spectrum 
with only a shallow structure, as shown by the dashed red line in Fig. [2^. The spectrum 
up to 1 TeV could be reasonably well described by a power line proportional to E in 
agreement with data from the H.E.S.S. experiment (High Energy Stereoscopic System, 
an earth-bound array of Cherenkov telescopes, which measures high energy gamma rays 
from the Cherenkov light in the air shower), as shown in Fig. [2jx Also these data show 
no indication of a structure in the electron spectrum, but rather a power-law spectrum 
with a spectral index of 3.0 ±0.1 (star.) ±03(syst.), which steepens at about 1 TeV flU, 
indicative of the cut-off of nearby astrophysical sources. 



CONSTRAINTS FROM ANTIPROTONS 

Charged particles frequently change directions in the irregular magnetic fields of the 
Galaxy. This kind of random walk is usually described by a diffusion equation. The 
diffusion equation describing the transport of cosmic rays (CRs) in the Galaxy is solved 
numerically in the publicly available GALPROP code [25 J. In its simplest version, 
sometimes called the conventional model, the diffusion coefficient is isotropic and has 
everywhere the same value and the spectral shape of the locally observed CRs is assumed 
to be representative for the Galaxy. Details about Galactic propagation models can be 
found in a review by Strong et al. E6l . 

The isotropic propagation models have a few important deficiencies: i) they do not 
allow for significant convection of CRs by gas streams, which are therefore assumed to 
have negligible velocities. However, this disagrees with the recent ROSAT data on X- 
ray spectra in the halo 11281 ii) they cannot explain the large bulge/disc ratio of the 511 
keV positron annihilation line, as observed by INTEGRAL, see Ref. [29 ] for a review 
iii) they tend to overproduce the number of gamma rays towards the Galactic centre 
in comparison to the gamma ray production away from the Galactic centre ll3~0"l . These 




FIGURE 3. Left: Antiproton/proton ratio in an isotropic propagation model from Donato et al. Il27l . 
Right: Antiproton spectrum in an anisotropic propagation model. 



deficiencies can be remedied if one allows for a different diffusion constant in the halo 
and the disc OTTl . In such a model the observed amount of antiprotons tends to be a factor 
of about two below the observed rate. Such a tendency was already observed before for 
the conventional GALPROP model ll32l and disagrees with semi-analytical calculations, 
which can describe the observed antiproton fluxes reasonably well, as shown in Fig. [3^. 
The antiproton production in an anisotropic propagation model is shown in Fig. [3]). The 
secondary antiproton production is roughly a factor two below the data, so one needs 
an additional source of antiprotons in such a model. This can be nicely provided by 
the annihilation of a 60 GeV neutralino, as shown by the red area. Note that low mass 
WIMPs have a spectral shape similar to the secondary (+tertiary) antiprotons, so their 
contributions cannot be distinguished from the spectral shape. In such models half of the 
antiprotons can be attributed to dark matter annihilation, thus relaxing the constraints on 
DMA models discussed by Donato et al. Il27l . 

INTERPRETATION 

The data described in the previous sections has caused an explosion of papers on possible 
interpretations in terms of new physics, varying from new astrophysical sources to dark 
matter annihilation. Unfortunately, none of the explanations can be excluded and what is 
worse, they all seem possible, so they can contribute simultaneously. For astrophysical 
sources pulsars are outstanding candidates for positrons, but also secondary acceleration 
of positrons in SNR may contribute or local nearby sources with correspondingly lower 
energy losses may harden the positron fraction. Also hadronic background may be non- 
negligible, so the total contribution to the positron fraction may be described by 

\Data >= abg\Backgr > +a p i s \Pulsar > +a mr \SNR > +ai 0C \Src > +a<i ma \DMA > 

(1) 
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FIGURE 4. CR spectra. Adapted from Lipari E). 



Of course the sum of the contributions cannot be more than the data, so the coefficients a 
in front of the contributions are constrained by unitarity. However, at present it looks like 
most contribution may saturate the data or not contribute at all, so most of the coefficients 
in Eq.[T]are allowed to have values between and 1. Clearly, more independent data are 
needed to constrain the different contributions. We will now discuss in more detail the 
various contributions. 



Background 

Separating antimatter from matter particles in cosmic rays (CRs) is not an easy task, 
because of the preponderance of protons, followed by electrons, which have already a 
flux, which is several orders of magnitude lower, as demonstrated in Fig. [4] [33l . 

To separate positrons from protons one needs a rejection of at least 1 : 10 5 , i.e. only 
one proton in 10 5 may simulate a positron, if one wants a background at the percent 
level. This can be achieved by a combination of a magnetic spectrometer measuring 
the charge, a calorimeter for measuring the shower profile and a transition radiation 
detector for measuring the velocity. This is the classical design used in the HEAT 
experiment and foreseen for the upcoming AMS-02 experiment. For PAMELA the 
same configuration was planned, but the TRD detector had to be abandoned in the 
last moment for technical reasons. Nevertheless, the fine granularity Si-W calorimeter 
provides typically a rejection of 10 -4 at 10 GeV to 10~ 5 at 100 GeV, as shown by the 
PAMELA beam test Poll . 

The only doubt coming to mind about the amount of background arises from the 




FIGURE 5. Left: AMS-02 hadronic background to the positron fraction for a proton rejection level 
of 1CT 5 . Note that the indicated DM model to describe the PAMELA data points has a similar slope 
as the hadronic background. From Battiston [34|. Right: the same conclusion is reached by the HEAT 
experiment: if the rejection is 3-10 4 , the positron fraction has the shape expected from the background. 
From Tarle |35|. 



fact, that for new physics the increase could have had any slope, but the increase is 
exactly what one would expect from an hadronic background, as shown by the estimates 
from AMS-02 and HEAT in Fig. [5J Therefore a check of the hadron rejection with an 
independent detector, like the TRD, would have been useful and one is looking forward 
to independent data by the AMS-02 experiment in 2010, which is scheduled to be flown 
to the ISS in the middle of 2010 according to the NASA program. 

The FERMI background rejection is an order of magnitude worse than for PAMELA, 
because it does not have a magnetic spectrometer to measure charge and momentum. 
The estimated rejection power is almost flat at a level of 10 3 below 200 GeV and 
improves to 10 -4 at 1 TeV (2]. The hadron rejection therefore needs hard cuts and relies 
heavily on MC simulation, but even if there is a residual hadronic background, the wide 
structure observed by ATIC should have been observed. In addition, H.E.S.S. did not 
observe a structure (Fig. [2j)). So the structure seems improbable, but a confirmation of 
the hard electron spectrum from FERMI by the PAMELA data in the region of overlap 
would be nice. 



Pulsars 



Pulsars are rapidly spinning, magnetized neutron stars, which emit electromagnetic 
radiation along the magnetic poles, as can be observed as flashes of light on Earth, when 
the poles are directed towards us. They are also prime candidates for sources of primary 
electrons and positrons, since the strong rotating magnetic fields generate photons by 
synchrotron radiation, which in turn generate electrons and positrons by interaction 
with the magnetic field and these electrons and positrons generate again photons by 
synchrotron radiation. So such a chain reaction delivers lots of positrons and electrons, 
but no antiprotons. A review about pulsars and how they can explain the data was given 
by Profumo IfTTTl . The main question is: how many charged particles can escape from 
the strong magnetic fields. This is unknown, since especially young pulsars are still 




surrounded by their nebulae, which trap the electrons and positrons by the magnetic 
fields and energy losses for about 10 5 years. Therefore one has to consider pulsars older 
than this. Assuming the escape rate and/or the number of old and nearby pulsars to 
be free parameters such pulsars can nicely explain the PAMELA and FERMI data, as 
shown in Fig. [6] IT5l . If only a few nearby sources contribute, one could try to search 
for an anisotropy in the arrival direction of the positrons. However, even with only a few 
sources the diffusive scattering washes out anisotropics efficiently. 



Positron acceleration ins SNRs 

The increase in the positron fraction cannot be explained by current sources and 
propagation models, since one always expects a decreasing fraction, as shown by the 
GALPROP curve in Fig. [T^ and discussed more generally e.g. by Serpico lfT2l . It was 
suggested by Blasi IfPfll that an additional source of positrons could be the secondary 
positrons generated and accelerated in old SNRs, especially if they are located inside 
the high density of molecular clouds, the birthplaces of stars. If located at a distance 
of 1-2 kpc, this mechanism would in addition explain the cut-off at about 1 TeV in 
agreement with the HESS data fl4), since sources further away would not have such a 
hard spectrum due to energy losses. This physical idea can nicely explain the data on 
positrons and electrons, as shown in Fig. [7J although the strength of the effect depends 
heavily on the evolution of SNRs, which is not well known. Of course, explaining 
the rise in the positron fraction by acceleration of secondaries would also increase the 
antiproton/proton fraction, in disagreement with the data shown in Fig.[TJ). However, the 
increase could happen above 100 GeV, as shown by the simple modeling of Ref. Il37l . 
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FIGURE 7. Left: electron/positron spectra assuming secondary acceleration in SNR. From Ref. [14]. 



Local Sources 

As mentioned above the electron spectra from local sources might be different from 
sources further away due to energy losses. It is clear that in the disc the sources are not 
homogeneously distributed, since most of them are in spiral arms. The nearest arm is 
about 1-2 kpc from the solar neighborhood, so this high density of sources relative far 
away will deliver a large fraction of electrons with a relatively soft spectrum. Sources in 
the disc and local sources will deliver increasingly harder spectra and the superposition 
of all three components can easily describe the FERMI electron spectrum, as shown in 



Fig. 10 1 [19]. Note that the positrons are produced by the CRs in the gas, so they come 
from all distances and have a smooth distribution with an intermediately hard spectrum. 
The resulting positron fraction in this model first increases and then decreases, similarly 
to what would be expected from dark matter annihilation, but different from most other 
astrophysical explanations. Clearly, data above 100 GeV will be very interesting to see 
if it decreases. 



Dark Matter Annihilation 

Dark Matter, if a thermal relic, must have annihilated or decayed in order to reduce 
the high number density in the early universe, where all particles had similar number 
densities. Nowadays the number density of baryons is 10 orders of magnitude below 
the photon number density, as measured by the temperature of the cosmic microwave 
background. Given that the DM density, given by the WIMP mass times the number 
density, is at most a factor 6 above the baryonic density, the number density of WIMPs 
must also be small compared with the photon density, hence the WIMPS must have 
annihilated or decayed, if their number density was similar to the photon density in the 
early universe. 

The lightest neutralino in Supersymmetry is a perfect WIMP candidate, not only 
because it is expected to be stable because of R-parity, but also its annihilation cross 
section is exactly in the right ballpark needed for the 23% leftover of DM energy in 
the universe. It annihilated preferentially into heavy quarks, like b-quarks, because the 
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FIGURE 8. The thermally averaged neutralino annihilation cross section < <7v > in the m§m\i<i plane 
for bb and ft final states. 
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FIGURE 9. As Fig. [§] but for it, tt and W+W~) final states, 
annihilation into lighter fermions is helicity suppressed, as can be seen from Figs.[8]and 

m 

The quarks fragment into hadrons with a well-known spectrum and fraction of an- 
tiprotons, positrons, neutrinos and photons, as measured e.g. at the LEP accelerator [38J, 
so if one observes a rise in the positron fraction, one expects also an increase in antipro- 
ton fluxes. Since the observed rise in the PAMELA positron fraction is not accompanied 
by a strong rise in antiprotons, neutralino annihilation is not a candidate for this rise, 
at least if one assumes a thermal history of the universe. In non-thermal universes dark 
matter can arise from the decay of heavy particles, such as the moduli present in string 
theories. In this case freeze-out occurs at much higher temperatures, i.e. at higher densi- 
ties, so one needs a higher annihilation cross section to get the low relic density observed 
today ll39l . In this case one does not need a large boost factor. The problem with the an- 
tiproton problem stays, so one has to invoke propagation uncertainties or assume that 
other sources contribute to the positron flux as well, as discussed before. 
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FIGURE 10. Left: Fit of astrophysical sources to the FERMI electron spectrum and the PAMELA 
positron fraction. From Ref. [19|. Right: Fit of a leptophilic DM candidate to the FERMI electron 
spectrum and the PAMELA positron fraction. From Ref. J9)- 



If one ignores astrophysical contributions and investigates the possibility that DMA 
is completely responsible for the leptonic excesses then one has to resort to leptophilic 
WIMP candidates. This can happen if the WIMPS decay into new states, which are light 
enough to prevent antiprotons in the final states, i.e. typical WIMP masses are below 
1 GeV. For example, Nomura and Thaler [HI propose DM fermions decaying into an 
axion and a scalar with the latter decaying again to axions. For an axion mass in the 
range 360 - 800 MeV antiproton production is forbidden. Alternatively, Arkani-Hamed 
et al. [7] propose a new force carrier of a few GeV, e.g. a pseudo-Goldstone boson, 
into which a heavy WIMP can annihilate. Resumming the ladder diagrams of the <j) 
exchange can enhance the annihilation cross section for non-relativistic particles, i.e. 
WIMPs in the TeV range, by orders of magnitude [|40l . thus providing an explanation 
for the large cross section needed to explain the PAMELA and FERMI excesses via this 
so-called Sommerfeld enhancement [41J. Such models could also help explaining the 
large bulge/disc ratio of the INTEGRAL positron annihilation line ||29l . However, as in 
the case for the PAMELA and FERMI excesses, a very natural astrophysical explanation 
exists: Galactic winds, as observed by ROSAT, convect the positrons from the disc 
into the halo, thus suppressing the positron annihilation in the disc OTTl . Nevertheless, 
independent of astrophysical explanations, models with new force carriers are able to 
explain nicely the excesses in positrons and electrons as well, as shown in Fig.flOb fl9). 



CONCLUSION 



The observed excesses in electrons and positrons by PAMELA, FERMI, HESS and 
ATIC have led to an interesting mix of proposals, varying from new astrophysical 
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FIGURE 11. Left: preliminary FERMI data on diffuse gamma rays compared with EGRET data in 
the mid-latitude range (10° < b < 20°) [5 |. Right: published FERMI data from VELA pulsar compared 
with EGRET data flfe). Clearly, the published FERMI data agrees better with the EGRET data than the 
preliminary diffuse gamma ray data in the region around a few GeV. 



sources to DM annihilation. Unfortunately, all seem possible, which indicates how 
difficult it is to extract conclusive evidence about the origin of these charged CRs. This 
difficulty is intimately related to the fact that charged particles do not point to their origin 
but change direction by the Galactic magnetic fields. Gamma rays do point to the origin, 
so point sources can be subtracted and the isotropic extragalactic contribution can be 
estimated as well. 

FERMI will hopefully soon publish their data on diffuse gamma rays with energies 
ranging from tens of MeV to hundreds of GeV, i.e. the range of interest for typical WIMP 
scenarios; preliminary data indicate that the GeV excess in the diffuse gamma rays seen 
by EGRET is not observed anymore (5), although the calibration via the VELA pulsar 
did not show a large discrepancy between EGRET and FERMI in the region of the GeV 
excess [|42l|. as shown in Fig. [TT] Of course, the background for the diffuse gamma ray 
data is larger than for the strong VELA pulsar, which may indicate that the background 
was not under control in the EGRET experiment, since with the higher spatial resolution 
of the FERMI tracker and the better MC simulation the background should be better 
under control in the FERMI data. 

Clearly more independent data are welcome, which hopefully will be provided with 
the AMS-02 spectrometer to be flown to the ISS in the middle of 2010. But above all it 
is hoped that the LHC will provide direct clues to the nature of DM by finding WIMP 
candidates, which can then be compared to the indirect and direct searches for DM. 

Acknowledgments. I wish to thank Pran Nath for organizing an excellent SUSY09 
conference and many colleagues for illuminating discussions. This work was supported 
by the Deutsches Zentrum fur Luft- und Raumfahrt (DLR). 



REFERENCES 



1. O. Adriani, et al., Nature 458, 607-609 (2009), [MB .4 995 



2. A. A. Abdo, et al., Phys. Rev. Lett. 102, 181101 (2009), |0 905 . 0025| 

3. J. Chang, et al., Nature 456, 362-365 (2008). 

4. F. Aharonian, et al. (2009), 0905 . 01051 

5. T. A. Porter (2009), |0907.0294| 

6. O. Adriani, et al., Phys. Rev. Lett. 102, 051 101 (2009), [0810 . 4 994| 

7. N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, Phys. Rev. D79, 015014 (2009), 

10810 . 07131 



Y. Nomura, and J. Thaler, Phys. Rev. D79, 075008 (2009), [0810 . 5397 



9. L. Bergstrom, J. Edsjo, and G. Zaharijas, Phys. Rev. Lett. 103, 031 103 (2009), (0905 . 0333| 

10. A. M. Atoyan, F. A Aharonian, R . J. Tuffs, and H. J. Volk, Astr.& Astroph. 294, L41 (2005). 

11. S . Profumo (2008), |08 12 . 4457] 

12. P. D. Serpico, Phys. Rev. D79, 021302 (2009), [08T~0 .4846 



13. D. Hooper, P. Blasi, and P. D. Serpico, JCAP 0901, 025 (2009), 0810 . 1527 



14. P. Blasi, Phys. Rev. Lett. 103, 051104 (2009), |0 903 . 2794[ 

15. D. Grasso, et al., Astropart. Phys. 32, 140-151 (2009 )TT0905 . 63 6| 

16. D. Chowdhury, C. J. Jog, and S. K. Vempati (2009), |0 90 9 . 1182[ 

17. H. Yuksel, M. D. Kistler, and T. Stanev, Phys. R ev. Lett. 103, 51101 (2009), [0810 . 27 84 

18. Y. Fujita, et al., Phys. Rev. D80, 063003 (2009), 0903.5298] 

19. N. J. Shaviv, E. Naka r, and T Piran (2009), |0 902 . 037 6[ 

20. L. Bergstrom (200 9),|0903 . 4849| 

21. D. Hooper (2009), |0901.4090| 

22. S. W. Barwick, et al., Astrophys. J. 482, L191-L194 (1997) ,|astro-ph/97 03192| 

23. M. Aguilar, et al., Phys. Lett. B646, 145-154 (2007), [ast r o - ph / 7 3 1 5 4 1 

24. D. Hooper, and S. Profumo, Phys. Rept. 453, 29-115 (2007), |hep-ph/0701197[ 

25. A. Strong, and I. Moskalenko, http://galprop.stanford.edu/manuals/manual.pdf (2006). 

26. A. W. Strong, I. V. Moskalenko, and V. S. Ptuskin, Ann. Rev. Nucl. Part. Sci. 57, 285-327 (2007), 
|astro-ph/0701517| 

27. F. Donato, D Maurin, P. Brun, T. Delahaye, and P. Salati, Phys. Rev. Lett. 102, 071301 (2009), 

10810 . 52911 

28. D. Breitschwerdt, Nature 452, 826 (2008). 

29. N. Prantzos, New Astron. Rev. 52, 457-459 (2008), [0809 . 2 4 91| 

30. D. Br eitschwer dt, V. A . Dogiel, and H. J. Volk, Astron. Astrophys. 385, 216-238 (2002), 

|astro-ph/0201345l 

31. I. Gebauer, and W. de Boer (2009), 0910.2027. 

32. I. V. Moskalen ko, A. W. Strong, and O. Reimer, Astron. Astrophys. 338, L75-L78 (1998), 

|astro-ph/ 9808084] 



33. P. Lipari, Pamela Physics Workshop, Rom (2009), |http : / /pamela . roma2 . inf n . it/ 
workshop09/ slides_WS2009/Lipari . pdf . 

34. R. Battiston, private communication (2009). 

35. G. Tarle, Talk at PPC09, Oklahoma City (2009). 

36. M. Boezio, et al., Astropart. Phys. 26, 111-118 (2006). 



37. P. Blasi, and P. D. Serpico, Phys. Rev. Lett. 103, 081103 (2009), [0904 .0871 

38. C. Amsler, et al., Phys. Lett. B667, 1 ( 2008). 

39. G. Kane, R. Lu, and S. Watson (2009), |0 906 . 4765] 

40. J. Hisano, S. Matsumoto, M. M. Nojiri, and O. Saito, Phys. Rev. D71, 063528 (2005), [hep^ph? 
104124031 

41. R. Iengo, JHEP 05, 024 (2009), |0902.0688[ 

42. A. A. Abdo, et al., Astrophys. J. 696, '1084-1093 (2009), 0812 .2960] 



